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1.  INTRODUCTION 


High-voltage  charging  is  now  recognized  (Parks  and  Katz,  1981;  Katz 
and  Parks,  1933)  to  be  an  operating  hazard  for  larger  spacecraft, 
including  the  Shuttle  crbiter,  passing  through  the  auroral  plasma  in  low 
polar  orbit.  It  is  important  to  develop  methods  for  predicting  which 
combinations  of  environmental  conditions  and  spacecraft  properties  will 
result  in  high-voltage  charging. 

The  work  presented  in  this  Report  is  in  two  parts.  Section  2 
contains  a  sim.ple  approximate  theoretical  prediction  of  the  required 
conditions  for  high-voltage  polar-orbit  charging.  The  results  of  this 
derivation  suggest  that  spacecraft  potentials  are  likely  to  depend  more 
strongly  on  the  ratio  of  am.bient  flux  of  high-energy  electrons  to  that  of 
all  ions,  than  on  any  other  environmental  parameter.  In  Sec.  3, 
calculations  are  made  of  secondary-electron  escape  currents  from  space¬ 
craft  surfaces,  as  influenced  by  magnetic  fields  having  various 
directions  relative  to  these  surfaces.  For  realistic  values  of  electron- 
repelling  surface  electric  fields,  the  results  show  an  extreme  sensitivity 
of  escaping  currents  to  small  changes  in  surface  orientation,  for 
surfaces  alm.cst  parallel  to  the  m.agnetic  field  direction.  This  implies 
that  the  occurrence  of  high-voltage  charging  in  marginal  circum.stances 
may  dsoend  ver^y  strongly  on  the  precise  orientation  of  a  surface. 
Apoendix  A  contains  a  listing  of  the  computer  program  used  to  perform 
these  calculations. 
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2.  ESTIMATE  OF  REQUIRED  ENVIRONMENTAL  CONDITIONS  FOR 
LOW-POLAR-ORBIT  CHARGING 


I 


m 


Ip.  this  Sectiop,  we  show  that  spacecraft  surface  potentials  are  likely 
to  depend  more  strongly  on  the  ratio  of  ambient  flux  of  high-energy 
electrons  to  that  of  ah  ions  than  on  any  other  applicable  environmental 
parameter.  To  do  this,  we  make  the  following  approximations. 

(1)  We  assume  that  magnetic-field  e.^fects  on  charged-particle  motion 
are  negligible.  This  assumption  should  be  acceptable  for  initial 
estimates  because  the  gyroradii  of  ions  and  high-energy  electrons  are 
generally  a  few  metres  or  larger,  especially  in  a  high-voltage  sheath 
(Laframboise,  1983,  Table  1),  and  collection  of  "cold"  (.vO.l  eV) 
ionospheric  electrons  by  a  negatively-charged  spacecraft  will  be  very 
sm.all,  so  their  density  is  well-approximated  by  a  Boltzmann  factor, 
independently  of  the  presence  of  a  m.agnetic  field. 

(2)  We  assume  that  ambient  high-energy  electrons  have  an  isotropic 
velocity  distribution.  Large  departures  from  this  have  been  observed 
in  auroral-plasma  conditions  (W.J.  Burke,  1984,  private  communic¬ 
ation),  but  this  should  not  seriously  affect  the  type  of  rough  estimate 
made  here.  Parks  and  Katz,  1981,  and  Katz  and  Parks,  1983, 
assum.ed  both  the  ion  and  electron  fluxes  to  be  unidirectional;  we 
discuss  this  ooint  later  in  this  Section. 


(3)  We  ignore  secondary-electron  emission;  magnetic-field  effects 
would  tend  to  suppress  this  on  some  parts  of  the  spacecraft  in  any 
case  (J.G.  Laframboise,  1983,  1985;  Sec.  3). 

(4)  We  assume  that  the  spacecraft  is  a  unipotential  sphere,  large 
compared  to  the  typical  ambient  Debye  length  of  5  1  cm.  We 
consider  only  overall  charging  of  the  spacecraft.  This  neglects  the 
possibility  that  local  high-voltage  charging  may  occur,  especially  on 
surfaces  in  the  spacecraft  wake. 

(5)  We  assume  that  both  ions  and  electrons  have  double-Maxwellian 
velocity  distributions,  with  the  colder  component  in  either  case 
having  a  temperature  of  0.1  eV,  and  the  hotter  1  keV  or  larger.  In 
the  spacecraft  reference  frame,  these  are  superposed  on  a  dri.^t 
velocity  equal  and  opposite  to  the  spacecraft  velocity. 

(6)  Ions  are  assumed  to  be  either  H  '  or  0  ' . 

Note  that  assum.ption  (3)  could  cause  a  false  prediction  that  high- 
voltage  charging  occurs,  while  assumption  (4)  could  cause  a  false 
prediction  that  it  does  not.  The  effects  of  assumptions  (1),  (2),  and  (5) 
are  less  clear;  these  could  conceivably  either  increase  or  decrease 
predicted  surface  potentials.  With  regard  to  (6),  assuming  that  the  ions 
are  H"*”  results  in  m.aximum  wake-filling  by  ions.  If  there  are  any 
electrically-isolated  surfaces  in  the  spacecraft  wake,  this  would  result 
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in  decreased  surface  potentials  (magnitudes);  assuming  O"*"  gives  the 
reverse. 

Probably  the  most  serious  difficulty  :n  formulating  a  theory  for 
Icw-crbit-charging  is  the  prediction  of  ion  collection  on  downstream 
surfaces.  As  mentioned  in  assumption  (3)  above,  we  avoid  this 
difficulty  by  considering  only  total,  rather  than  local,  ion  collection,  on 
a  unipotential  sphere.  Kanal  [1962,  Eq.  (63)]  gives  an  expression  for 
the  icn  current  collected  by  such  a  sphere  from  a  drifting  Maxwellian 
plasma  in  the  lim.it  of  zero  potentials  (relative  to  space  potential),  as 
follows: 

i,  =  (S.  +  ^)  erf(Sp  +  exp(-S.‘)]  (2.1) 

where  i.  =  h/L;,  L;  is  the  ion  random  current  en.  (kT;/27rm.)^, 

1  I  01  01  100  I  X 

S.  =  U/(21<T./m;)^  is  the  ion  speed  ratio,  U  is  the  ion  drift  speed 
relative  to  the  spacecraft,  e  is  the  magnitude  of  the  electronic  charge,  k 
is  Boltzmann’s  constant,  and  m^  T.,  and  n.^  are  ion  m.ass,  temperature, 
and  am.bient  number  density.  We  assum.e  that  U  =  8  km./sec,  corres¬ 
ponding  to  low  circular  orbit. 

We  need  to  take  account  of  the  effect  of  a  large  ion-attracting  surface 
potential  on  ion  collection,  in  the  limit  of  sm.all  Debye  length  Xq 
com, pared  to  the  sphere  radius  r  .  To  do  this,  we  use  a  result  of 
Parrct  et  al  (1982).  These  authors  show  that  for  a  probe  in  a  collision¬ 
less,  nonm.agnetized,  Vlaxwellian  plasma  having  =  1  and  without 


Ion  drift,  and  in  the  limit  when  0  but  -eOg/kT  >>  i  [where 

these  limits  must  be  approached  in  such  a  way  that  (-e0^/kT)  (X^/r^) 

remains  <<"  1,  i.e.  sheath  thickness  remains  <<  sphere  radius],  the  ion 
(attracted-particle)  current  is  larger  than  the  random  current  by  a  factor 
of  1.45.  This  factor  represents  the  effect  of  the  "presheath"  potential  on 
ion  collection.  Even  though  several  of  their  assumptions  are  unfulfilled 
in  our  case,  the  resulting  effects  on  ion  collection  are  probably  small 
enough  for  our  purposes.  We  therefore  multiply  Eq.  (2.1)  by  the  same 
factor  to  obtain  an  estimate  of  total  ion  collection  as  influenced  by 
surface-potential  effects.  The  resulting  ion-current  dependence  on  ion 
speed  ratio  is  plotted  in  Fig.  2.1.  For  O'  ions  at  T  =  0.1  eV  (1160K), 
H  at  0.1  eV,  O '  at  1  keV,  and  H '  at  1  keV,  S.  =  7.31,  1.83,  0.0731, 
and  0.0183  (the  latter  two  are  effectively  zero),  respectively.  The 
corresponding  ion-current  enhancement  factors  (values  of  i^  from  Fig. 
2.1  a.-^e  9.50,  2.59,  1.45,  and  1.45,  respectively. 


If  the  ambient  ions  are 


,  the  ion  collected  current  is  now  given  by; 


I.  =  47rr  ^  en. 

1  s  :c 


2nm. 
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(2.69) 
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en 


ih 


kT. 


ih 


:7rm. 


(1.45) 


(2.2) 


where  the  subscripts  ic  and  ih  refer  to  the  cold  and  hot  ion  populations. 
If  the  ions  are  O'*’,  then  the  factor  2.69  in  (2.2)  should  be  replaced  by 
9.50. 
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.he  electron  collected  ccrrent  is: 


,  ,  I* 


(2.3) 


-  -eh  t  j 


1:  high-voltage  charging  occurs,  then  ~^<P^  >>  '^‘c.n‘  first  term  cn 

-.ha  right-hand  side  of  this  equation  becomes  negligible. 


rcr  current  balance,  1.  =  1  .  ;his  leads  to: 
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vhere  vm../m  =  *43  for  H'  ions,  .herefore: 
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•  H  ions,  with  43  and  2.69  replaced  by  172  and  9.50  for  O'  ions. 
:s  is  ecuivalent  to: 


hot-electron  am.bient  flux 


.  /'t^T  ,  =  jn  i  2.69  (coid-ion  ambiera  fluxy  -r  i.45  (riot  ion 
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i,3h-voit:age  crargir.g  to  become  probable,  the  argum.ent  of  the  in 
m.'jGt  be  close  to  or  larger  than  e  111  2.72,  i.e: 


_ l-Qt-aleotrcr.  am.bient  flux _ 

2.55  'oold-icn  am.bient  flux)  +  i  .'15 (hct-icn  ambient  flux) 


2.72. 


.-or  O  /H  mixtures  and  for  hot-ion  temoeratures  other  than  1  keV. 


:  this  result  is  straightforward.  Since  any  hot  ions  are 


. iKe.v  to  rave  .../.. 

•  n  • : 


iC'*.  the  hot-ion  ambient  flux  will  exceed  the 


old-:or  ambient  flux  i:  the  net  ions  constitute  m.ore  than  about  1%  of  the 
ttal  am.bien.t-ion.  "um.ber  density.  Equation  (2.7)  indicates  that  the 


tnset  of  h:gh-voltag9  charging  can  be  expected  to  depend  primarily  on  the 
ra':o  or  hot-electron,  am.bient  flux  to  the  ambient  flux  of  all_  ions,  as 
m.en.tioned  at  the  begm.nin.g  of  this  Section.  This  com.pietes  our  argument 

a  .  ■  w  i  ■y  I  .  «  W  W  I  <  I  k  • 


In  analvzing  spacecraft  data,  one  is  therefore  likely  to  find  better 
torrelation  cf  spacecraft  voltages  with  the  ratio  which  appears  in  Eq. 
’2."''  ,c.r  som.e'hin.E  nearly  equal  to  it)  than  with  any  other  m.easurable 
tuentity,  such,  as  electron  cr  ion  density  or  average  energy,  taken 
.-div'.dually.  In  calculating  values  of  this  ratio,  the  amfoient  fluxes 
v'-ich  are  involved  need  to  have  been  measured  simultaneously  on  the 
tam.e  scacecroft.  Even  though  the  apprcxim.ations  m.ade  in  deriving  (2.7) 
ar^  severe,  and  the  precise  dependence  cf  spacecraft  voltages  on  this 
~atio  m.ay  therefore  differ  substantially  from  that  given  in  Eq.  (2.7) 
a-cl/or  the  coefficients  in  the  ratio  itself  m.av  need  to  be  modified),  our 


.M'.'b-'-r.'-N' 
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general  ccnclusion,  i.e.  that  spacecraft  voltages  should  correlate  most 
strongly  with  this  ratio  (or  something  nearly  eoual  to  it),  is  likely  to 
rem.ain  va'id.  Furtherm,ore,  the  dependence  of  spacecraft  voltages  on 
this  flux  ratio  is  likely  to  retain  an  approxim.ately  exponenti"'!  form,.  In 
situations  where  m.cst  secondary  and  backscattered  electrons  emitted  by 
the  spacecraft  will  escape  (see  Sec.  3),  prim.ary-electron  incident  fluxes 
will  be  approxim.ately  cancelled  for  .m.any  spacecraft  m.aterials  by 
electron  escape  at  incident  energies  up  to  a  .^'ew  keV  (La.^ramboise  et  al, 
1982;  La.^ram.boise  and  Kamitsum.a,  1983),  so  the  hct-electrcn  am.bient 
flux  term  :n  (2.7)  .needs  to  be  modified  accordingly. 

The  most  serious  approximation  .m.ade  in  deriving  (2.7)  is  probably 
ite.m.  (4)  in  the  list  at  the  beginning  of  this  Section.  This  is  because  ion 
fluxes  on  downstream,  surfaces  are  likely  to  be  very  m.uch  smaller  than 
their  average  over  the  entire  spacecraft.  They  are  also  likely  to  be 
strongly  dependent  on.  spacecraft  geometry,  local  surface  potential 

^  J- 

distribution,  and  O  /H  concentration  ratio.  Therefore,  the  critical 
value  of  ambient  flux  ratio,  at  which  the  onset  of  high-voltage  charging 
occurs,  is  likely  to  vary  substantially  among  spacecraft  having  different 
geometries  and  surface  .materials.  In  particular,  for  spacecraft  having 
electrically-isclated  downstream  surfaces,  this  critical  ratio  is  likely, 
because  of  local  changing  on  these  suri'aces,  to  be  much  lower  than  for 
spacecraft  which  have  an  entirely  conductive  surface. 

Funtherm.cre,  in  contrast  with  the  situation  for  total  ion  collection, 

8 


■  ■' ‘S'.  ' ’ 

W*" W V wVm..  y./ m-rv  mJ'Zi 


there  is  no  known,  simple,  reliable  method  for  estimating  ion  fluxes  on 
downstream  surfaces.  Parks  and  Katz  (1983a, b)  have  developed  an  ion 
flux  calculation  for  the  downstream  point  on  a  sphere  in  a  potential 
which  has  a  given,  simple  analytic  form.  Detailed  numerical 
simulation,  which  includes  realistic  self-consistent  spacecraft  sheath 
potential  distributions,  and  which  probably  needs  to  involve  at  least 
some  ion  orbit-following,  therefore  appears  to  be  essential. 

Preliminary  indications,  from  work  of  this  type  presently  in  progress 
(L.W.  Parker  and  J.G.  Laframboise,  to  be  published),  are  that  conditions 
on  the  "shoulder"  or  "side-point"  regions  of  spacecraft  (surface  material 
and  geometry  of  this  region;  local  surface  curvature  is  probably 
important)  may  strongly  influence  potentials  of  downstream  surfaces, 
because  of  detailed  effects  on  ion  trajectories.  The  geometry  being 
investigated  first  is  an  infinite  nonconducting  cylinder  transverse  to  the 
ion  drift  direction,  and  preliminary  results  also  indicate  that  the 
location  of  m.axim.um  negative  surface  voltage  is  not  at  the  downstream 
point  but  at  two  points  symmetrically  located  on  either  side  of  it. 

So  far,  we  have  not  mentioned  the  difficulties  which  can  arise  in 
measuring  the  ambient  ion  fluxes  whicfi  appear  in  Eq.  (2.7).  So  far,  we 
have  also  defined  "ambient  flux"  to  be  that  measured  in  an  Earth-fixed 
reference  fram.e.  The  alternative  would  be  to  define  it  as  that  measured 
in  the  spacecraft  frame,  i.e.  including  ram  effects.  Ion  fluxes 
measured  by  spacecraft  instruments  are  strongly  influenced  by  ram 
effects.  In  fact,  the  numerical  factors  2.69,  1.45,  and  9.50,  which 


appear  in  Eq.  (2.7)  and  the  associated  discussion,  already  constitute  a 
rough  ram-effect  correction,  but  for  total  current  to  a  sphere,  not  for 
local  collection  by  a  forward-facing  instrument  apertu.''e.  It- may  happen 
that  the  ram-effect  correction  factors  for  an  instrument  are  nearly 
equal  to  the  above  factors,  so  that  the  instrumient  measurement,  without 
any  correction,  already  gives  a  good  estimate  of  the  denom.inator  of 
Eq.  (2.7).  In  any  case,  the  response  of  the  instrument  will  depend  on  its 
geom.etry,  and  this  problem  has  already  been  treated  by  other  authors 
(Parker,  197Q;  Parker  and  Whipple,  1970;  Whipple  et  al,  1974;  Chang 
et  al,  1979;  Singh  and  Eaugher,  1981;  Com.fort  et  ai,  1982; 
Laframboise,  1983),  so  we  do  not  discuss  it  here. 

Parks  and  Katz  (1981)  and  Katz  and  Parks  (1983)  have  estimated 
charging  potentials  on  spherical  spacecraft  of  0.5m  and  5m  radius, 
assuming  that  the  ions  are  O^,  the  hot  electron  temperature  is  5 
keV,  and  spacecraft  speed  is  8  km/sec.  Their  results  can  be  compared 
directly  with  those  given  by  our  Eqs.  (2.5)  -  (2.7).  They  have  used  the 
theory  of  Langmuir  and  Blodgett  (1924)  to  obtain  values  for  sheath 
radius  as  a  function  of  spacecraft  potential.  They  present  spacecraft 
potentials  as  functions  of  the  ratio  k  of  hot  ("precipitating")  electron  ram 
current  to  ion  ram  current.  To  make  a  comparison,  their  value  of  k 
needs  to  be  expressed  in  term.s  of  our  ambient  flux  ratio.  They  have 
assumed  the  ambient  electron  flux  to  be  unidirectional.  To  convert  to  an 
equivalent  isotropic  flux,  we  note  that  current  to  a  sphere  =  ^ 

isotropic  (random)  flux,  but  =  x  unidirectional  (ram)  flux. 
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n 

•* 

n 

hererore.  equivalent  isotropic  rlux  -■  i  x  unidirectional  flux,  for  a 

Also  for  a  sphere,  the  ratio  of  ion  ram  to  random  currents  is 

j/iSkT./ irm  j ^  ;  Vtt  S..  Using  S.  =  7.31,  this  ratio  =  6.48,  so 

herefcre; 

^  _  hot  electron  ram  current 
"  6.43  X  total  ion  random  current 

rrr  ^  X  hot  electron  ram  flux 


5.48  X  47rr  ^  x  total  ion  random  flux 


1_  ^  hot  electron  (equivalent)  random  flux 
5.43  total  ion  random  flux 


.  w; 


6740 


X  our  flux  ratio  R. 


Vith  coefficients  for  O'  used,  our  Eq.  (2.6)  gives: 


o„  =  -5000  In  (R/9.50). 


(2.9) 


hgure  2.2  shows  our  result  and  theirs  [from  their  Fig.  3  (1981)  or 
hg.  2  (1983)],  plotted  together.  At  larger  potentials,  the  combined  set 
f  results  shows  a  monotonic  progression  toward  increased  charging  for 
arger  spacecraft.  For  -cp^  <  350V,  their  5m  sphere  shows  more 
hanging  than  our  large-radius-limit  sphere.  This  is  because  their  ion- 
u'rrent  enhancement  factor,  which  is  determined  by  the  size  of  a  sharp- 
cged  L.angmuir-Blodgett  sheath,  falls  below  ours,  which  includes  the 
ffect  of  a  quasineutral  presheath.  This  discussion  suggests  that  the 


3 


tendency  toward  high-voltage  charging  always  increases  with  spacecraft 
size,  but  magnetic-field  effects  may  modify  this  (Laframboise,  1983, 
Sec.  1).  The  corresponding  curves  for  local  charging,  on  surfaces  in  a 
spacecraft  wake,  will  lie  to  the  left  of  those  shown  in  Fig.  2.2,  but 
these  remain  to  be  computed  numerically. 


3.  CALCULATION  OF  SECONDARY-ELECTRON  ESCAPE  CURRENTS 
FROM  NEGATIVELY-CHARGED  SPACECRAFT  SURFACES  IN  A 


I 


'v^ 


MAGNETIC  FIELD 


3.1  SUMMARY 


In  low  Earth  orbit,  the  geomagnetic  field  is  strong  enough  that  secondary 
electrons  emitted  from  spacecraft  surfaces  have  an  average  gyroradius  much 

smaller  than  typical  dimensions  of  large  spacecraft.  This  implies  that  escape 
of  secondaries  will  be  strongly  inhibited  on  surfaces  which  are  nearly  parallel 
to  even  if  a  repelling  electric  field  exists  outside  them.  This  effect  is 
likely  to  m.ake  an  important  contribution  to  the  current  balance  and  hence  the 
equilibrium  potential  of  such  surfaces,  making  high-voltage  charging  of  them 
more  likely.  We  present  numerically-calculated  escaping  secondary-electron 
fluxes  for  these  conditions.  For  use  in  numerical  spacecraft-charging 
simulations,  we  also  present  an  analytic  curve-fit  to  these  results,  accurate  to 
within  3%  of  the  emitted  current. 


3.2  INTRODUCTION 


The  prediction  of  high-voltage  charging  or  other  environmental  effects  on  a 
spacecraft  in  low  Earth  orbit  appears  likely  to  be  more  complicated  than  in 
geostationary  orbit,  for  at  least  three  reasons. 


These  reasons  are:  (a)  space-charge  effects  (on  sheath  and  wake  potentials) 
are  more  important,  because  space-charge  densities  are  much  higher  (the 
Debye  length  is  no  longer  >>  typical  spacecraft  dim.ensions)  (b)  ion  flow 
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effects  are  more  important,  because  spacecraft  orbital  speed  >  ion  thermal 
speeds  (c)  the  geomagnetic  field  is  likely  to  have  an  important  influence  on 
charged-particle  motions  because  is  now  much  larger,  and  not  all  of  the 
average  particle  gyroradii  of  importance  are  any  longer  >>  typical  spacecraft 
dimensions. 

We  wish  to  investigate  an  important  consequence  of  (c),  which  concerns  the 
escape  of  secondary  electrons  emitted  from  spacecraft  surfaces.  Our  discus¬ 
sion  will  also  apply,  with  minor  modifications,  to  photoelectron  or 
backscattered-electron  escape.  In  low  Earth  orbit,  in  the  auroral-zone 
geomagnetic  field  (|B^j=0.44  gauss  =  4.4  x  10  T),  the  gyroradius  of  a 
"typical"  3eV  secondary  electron  and  a  10  keV  auroral  electron  are  13  cm  and 
8  m,  respectively.  The  average  gyroradius  of  "cold"  ionospheric  electrons 
(temperature  T  =  0.1  eV)  in  the  same~S  is  even  smaller  (2  cm),  but  this  is 
not  an  im.portant  param.eter  in  most  cases  because  these  electrons  are  repelled 
if  the  spacecraft  potential  is  negative,  and  their  density  is  then  well- 
approximated  by  a  Boltzmann  factor,  which  is  unaltered  by  B*”  effects. 

The  reason  why  B^ affects  secondary-electron  escape  is  shown  in  Fig.  3.1. 

In  Fig.  3.1(a),  the  spacecraft  surface  is  perpendicular  to  B^,  and  the  emitted 
electrons,  which  experience  an  electric  force  -e  directed  away  from 
the  surface,  all  escape,  helping  to  discharge  it.  In  Fig.  3.i(b),  the  spacecraft 
surface  is  nearly  parallel  to  B^,  and  almost  all  of  the  emitted  electrons  return 
to  it,  even  though  they  still  experience  an  electric  force  directed  away  from  it. 
These  electrons  therefore  are  unable  to  help  discharge  it,  so  a  surface  nearly 
parallel  to  B^  is  more  likely  to  charge  to  a  large  negative  voltage.  Note  that 
the  component  of  which  is  perpendicular  to  B^results  only  in  an  B^<  B^’drift 
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parallel  to  the  surface. 
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For  any  object  much  larger  than  13  cm,  the  escape  of  secondary  electrons 
will  be  strongly  affected  by  this  process.  For  example,  most  surfaces  on  the 
Shuttle  are  effectively  "infinite  planes"  by  this  criterion.  On  the  other  hand, 
the  average  gyroradius  of  high-energy  auroral  electrons  is  comparable  to 
Shuttle  dimensions,  so  the  deposition  of  these  electrons  onto  Shuttle  surfaces 
is  likely  to  be  only  moderately  inhibited. 


For  a  larger  object  (size  >>  8  m),  deposition  of  auroral  electrons  will 
also  become  strongly  orientation-dependent,  with  both  collection  and  escape  of 
electrons  now  being  inhibited  on  surfaces  nearly  parallel  to  This  suggests 
that  high-voltage  charging  of  such  surfaces  may  be  more  likely  on  objects  of 
intermediate  size  than  on  either  larger  or  smaller  ones.  In  the  calculation  of 
Parks  and  Katz  (1981),  Katz  and  Parks  (1983),  the  tendency  toward  high- 
voltage  charging  increased  with  spacecraft  size  because  in  their  model,  ion 
collection  increased  less  rapidly  with  spacecraft  size  than  did  electron 
collection.  To  determine  which  of  these  two  effects  predomiates  will  require 
more  detailed  calculations  than  have  been  done  so  far. 


As  already  mentioned,  strong  ion  flow  effects  also  are  generally  present  in 
low  orbit;  the  ion  speed  ratios  (flow  speed/most  probable  ion  thermal  speed) 
for  H  at  1  keV,  H'*'  at  0.1  eV,  and  at  0.1  eV  are  0.02,  1.8,  and  7.3, 
rescectively.  Whenever  the  latter  is  the  predominant  ion  species,  ion 


collection  on  downstream  surfaces  will  therefore  be  strongly  inhibited.  If  a 
surface  is  simultaneously  downstream  and  nearly  parallel  to  ^ as  is  likely  to 
be  the  case  in  the  auroral  zones,  then  the  tendency  for  high-voltage  charging  to 
occur  on  it  will  be  greatly  increased  (Fig.  3.2). 


To  "straightforwardly"  include  B*^effects  on  secondary-electron  emission  in  a 
large  two  or  three  dimensional  simulation  program  would  involve  the 

r.um.erical  integration  of  very  large  num.bers  of  secondary-electron  orbits.  The 
resulting  ccm.puting  costs  usually  would  be  formidable,  especially  since  these 
orbits  would  have  relatively  large  curvatures,  A  desirable  alternative  is  to 
"parameterize'  the  situation  by  treating  in  advance  a  simplified  but  still 
sufficiently  realistic  model  problem.  In  order  to  do  this,  we  make  the  approx¬ 
imations  described  in  Sec.  3.3. 

3.3  THEORY  FOR  E*’N0RMAL  TO  SURFACE 

We  assume  that  the  spacecraft  surface  is  an  infinite  plane,  and  the  electric 

and  magnetic  fields  and  B*"  outside  it  are  uniform.  In  the  work  presented 

here,  we  also  assumie  that  the  electric  force  -eB^  on  electrons  is  directed  along 
the  outward  ncrm.al  to  the  surface;  here  e  is  the  magnitude  of  the  elementary 

charge.  This  assumption  is  to  be  relaxed  later  (J.G.  Laframboise,  to  be 
published)  in  order  to  per.mit  variations  of  potential  along  the  surface  to  be 
taken  into  account.  We  assume  that  the  secondary  electrons  are  emitted  with 
a  Vlaxwelliar.  distribution  corresponding  to  a  temperature  T.  The  ratio  1=1/1^ 
of  escaping  to  emitted  flux  is  then  a  function  of  two  parameters:  the  angle  0 
between  the  surface  normal  and  the  direction  of  (Fig.  3.3),  and  a  parameter 
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describing  the  strength  of  E.  A  convenient  choice  for  this  parameter  is  the 
difference  in  potential  across  a  mean  secondary-electron  gyroradius  a  = 
(1/eB)  (TrmkT/2)  ^  ,  divided  by  kT/e,  where  m  is  electron  mass  and  k  is 
Boltzmann’s  constant. 

This  quotient  is: 


where  E  =  i  and  B  =  iB*"  I . 

This  quantity  also  has  an  alternative,  more  useful  interpretation:  it  is  the 

ratio  of  the  magnitude  |E^  of  the  drift  speed,  to  one-half  the 

mean  thermal  speed  (SkT/rrm)^  of  the  emitted  electrons.  It  is  useful  to 

estimate  the  value  of  e  for  a  high-voltage  spacecraft  sheath  in  low-orbit 

conditions.  To  do  this,  we  use  the  sheath  solution  of  Al’pert  et  al  (1965, 

Table  XXIV  and  Fig.  72).  For  a  1  kV  and  a  5  kV  sheath  around  a  sphere  of 

radius  3m.  in  a  collisionless  plasma  having  an  ambient  ion  temperature  of 

_3 

O.IV,  number  density  of  3  x  IQS  cm  ,  and  resultant  (ion)  Debye  length  of 
0.43  cm,  their  results  give,  respectively,  sheath  thicknesses  of  2.6  and  6.1  m, 

-S 

and  surface  electric  fields  E  =  0.86  and  2.9  kV/m.  Using  B  =  4.4xl0  T 
and  T  =  3  eV  fcr  secondary  electrons,  we  then  obtain  e  =  33.9  and  114.2. 

Both  of  these  are  relatively  large  values,  whose  significance  can  be  understood 
if  we  consider  what  would  happen  if  e  were  infinite. 


In  this  limit,  it  is  easy  to  show  that  secondary  electrons  would  all  escape 
unless  were  exactly  parallel  to  the  surface  (6  were  90°).  This  can  be 
shown  as  follows.  In  this  limit,  secondary  electrons  would  have  no  "thermal" 
m.otion.  The  (y,z)  projection  of  their  .motion  would  then  be  similar  to  that 
shown  in  Fig. 3. 4  .  This  motion  would  be  the  sum  of:  (i)  an  "Adrift  in  the 
y  direction  (ii)  a  uniform  acceleration  along  B  ,  whose  projection  in  the  (y,z) 
plane  would  be  upward  (iii)  just  enough  gy.^omotion  to  produce  a  cycloidal  path 
when  com.bined  with  (i),  so  that  in  the  absence  of  (ii),  the  electron  would  (just) 
return  to  the  surface  at  the  end  of  each  gyroperiod.  In  the  presence  of  (ii), 
these  "return  points"  are  displaced  upward  by  progressively  increasing  amounts 
(Fig. 3. 4),  so  the  electron  can  never  return  to  the  surface,  unless  is  exactly 
parallel  to  the  sur.'^ace,  so  that  the  upward  component  of  -eE^  along  vanishes. 
If  -eb^  has  a  component  parallel  to  the  surface,  this  conclusion  needs  to  be 
modified  (J.G.  Laframboise,  to  be  published). 

This  result  suggests  that  for  large  f^inite  values  of  e  (including  the  values 
calculated  above),  electron  escape  is  likely  to  be  almost  complete  except  for  0 
very  near  90°,  where  it  should  drop  to  zero  very  steeply.  The  occurrence  of 
high-vcltage  charging  in  marginal  circumstances  may  therefore  depend  very 
strongly  cn  the  precise  orientation  of  a  surface. 

The  escaping  secondary-electron  flux  is  given  by: 
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where:  is  the  initial  velocity  of  an  emitted  electron,  f(v^)  =  d  n/d  is 

the  velocity  distribution  of  emitted  electrons,  n  is  a  reference  number  density, 

and  H('^)  is  equal  to  1  for  escaping  electrons  and  0  for  those  which  return  to 

the  surface.  The  emitted  flux  is: 


=  n(kT/27rm) 


(3.3) 


We  also  introduce  the  dimensionless  velocity: 
lT  =  ^  (m/2kT)^  . 


Equation  (3.2)  then  becomes: 
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which  is  in  a  form  suitable  for  num.erical  summation.  The  quantities  u,.  <  , 

*  ilm,  i. 


j,.  .'j,.  .  are  the  vaiues  c:  :or  'vnicn  i-.  cnanpes  cetweer  U 

-im.Z  32  ° 

rrax 

=nci  i  for  eacr.  u  and  u  .  These  vaizes  .T.zst  be  found  bv  numerical Iv 
ox  oy  '  ' 

fatermining  which,  oarticle  orbits  reimcact  the  surface.  These  orbits  can, 
'.owever,  be  determined  in  analytic  form.,  wit.h  time  as  a  parameter.  To  co 
‘.his,  we  use  the  coordinate  system  shown  in  Fig.  3.3,  together  with  a  y-axis 
.not  shown)  directed  into  the  plane  of  the  Figure.  The  equation  of  motion  for 
an  electron  is: 


=  -  -  'zT*'-  X  3^ ^ . 


(3.6) 


We  solve  this  with  the  initial  conditions  ^  =  y  =  n  =  0,  Vj.  --  v  =  ,  a 

V  =  V  .  We  introduce  the  dimensionless  variables: 
on 


X  =  x/a.  V  =  v/a.  etc; 
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.n  the  present  v/ork,  e^  and  are  both  zero,  but  for  future  use,  we  "lave 
'etained  these  quantities  in  the  formulas  below.  We  obtain; 


u  ^  =  u  sin  6  +  u  cos  0; 
o(  ox  oz  ’ 


u  =  -u  cos  0  +  u  sin  0; 
or]  ox  oz  ’ 


f-  ~  ^  2  i  ^ 

<  =  — -t  —  u  t; 

y  =  Its  "=y  ■  ?  ^  J  ^  It 


'2  ^2  1  ,  ^  2 
yn  ^  y j  TT  q  ' 


r  2  ^21. 

n  =  —  '  f  sin  r 

‘  !_  yn  oq  ^  yj 
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by --iq}  I V’  f^-s; 


z  -  ^  cos  d  +  q  sin  0. 

Equations  (3.8)  can  also  be  differentiated  to  find  dz/dr.  The  numerical 
procedure  for  finding  the  quantities  ^  in  Eq.  (3.5)  then  involves 

'V  -V 

calculating  z  and  dz/dr  at  a  succession  of  points  along  an  orbit  (if  -eE  is 

normal  to  the  surface, the  electron  will  reimpact  during  the  first  gyropericd 

3  ''  0  ^  2n  if  at  all,  so  this  interval  always  suffices),  and  making  the 

appropriate  tests  on  these  quantities  to  find  out  whether  the  orbit  reim.pacts 

or  escaoes.  For  each  u  .and  u  .  this  is  done  for  a  succession  of  values 

cx,i  cy,j 

of  These  tests  also  yield  the  locai  minimum  of  z  (r)  if  one  exists. 

Wi'e^.ever  a  change  occurs  between  no  escape  and  escape  from  one  such  value 

0*  to  the  next,  an  interpolation  using  these  minima  can  be  used  to  provide 

‘ne  ccrresponding  '.oiue  of  u, .  ,  In  cases  where  they  are  unavailable,  the 

arithmetic  m.ean  of  the  two  succe^mive  u  values  is  used. 

oz 
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This  completes  the  definition  of  the  procedure  used  for  calculating  the  ratio 
l/'I^of  escaping  to  emitted  flux. 


3.4  RESULTS  AND  DISCUSSION 

Escaoing  secondary-electron  current  densities,  com.puted  as  described  in 


Sec.  3.3,  are  shown  in  lable  1  and  .^ig.  3.5.  Each  value  of  i  =  was 
calculated  using  191808  orbits,  evenly  spaced  in  the  intervals  -4.5  ^  u^^  ^ 


5,  -4.5  ^  u  si  4.5,  and  0  u  ^  4.5,  with  points  on  the  orbits  calculated 
’  oy  oz  ’  ^ 

intervals  Ar  =  7T'/45.  For  8  values  each  of  e  and  0,  the  resulting 


calculation  took  83  hr  total  on  a  Hewlett-Packard  lOOOF  minicomputer  v/ith 


Vector  Instruction  Set.  The  results  are  accurate  to  within  about  0.5%  or 


etter.  The  result  for  e  =  0  is  just  the  analytic  result  i  =  cos  0.  To  see  why 
this  is  so,  we  consider  the  electron  orbit  shown  in  Fig.  3.6,  which  has  been 
fictitiously  extended  so  as  to  pass  through  the  surface  and  re-emerge  from  it. 
In  the  absence  of  an  electric  field  (e  =  0),  this  orbit  has  the  same  speed  at  the 
re-em.ergence  point  C  as  at  the  emission  point  A.  Since  we  have  also  assumed 
that  the  em.itted  velocity  distribution  is  isotropic,  and  there.^ore  a  function  of 
speed  only,  the  real  orbit,  for  which  C  is  the  emission  point,  must  carry  the 
sam.e  population  as  would  the  fictitious  re-emerged  orbit.  The  flux  crossing 
the  reference  surface  DE,  which  is  i  is  therefore  the  same  as  if  such 


passages  and  re-emergences  actually  occurred,  and  is  the  same  as  if  another 
reference  surface  FG,  also  1  were  emitting  electrons  having  the  same 


•velocity  distribution.  However,  in  reality,  the  electrons  come  from  the  real 
surface  HJ,  which  is  not  i  and  all  the  electron-orbit  guiding  centers  which 
are  inside  any  given  m.agnetic-flux  tube  through  DE  will  also  be  inside  the 


projection  of  the  same  flux  tube  onto  HJ,  and  the  ratio  of  the  intersection 
areas  of  this  tube  with  HJ  and  DE  is  just  sec  0.  The  ratio  of  escaping  to 
emitted  flux  must  therefore  be  the  reciprocal  of  this,  or  cos  0,  as  stated 
above. 

Also  evident  in  Fig.  3.5  is  the  fact,  mentioned  in  Sec.  3.3,  that  when  e  is 
large  enough,  electron  escape  becomes  essentially  complete  except  when  0  is 
very  nearly  90°.  In  a  real  situation,  would  not  be  uniform,  but  would 
decrease  with,  distance  fro.m  the  surface,  contrary  to  our  assumptions.  Our 
results  can  therefore  be  expected  to  overestimate  electron  escape.  This  would 
probably  not  be  a  large  effect,  but  this  presum.pticn  rem.ains  to  be  verified.  An 
approximate  compensation  for  it  can  be  made  by  calculating  e  using  an  electric 
field  value  which  is  averaged  over  the  ."irst  mean  gyroradius  distance  from  the 
surface. 

The  results  in  Table  1  are  approximated  to  within  2.5%  of  by  the 
em.pirical  .‘^ormiula: 


a  =  1  +  1.35e‘-‘3°‘‘  exp^O. 083725  |l+tanh  [1.9732  in  (y^)] 


-0.07825  in  [1  +  (e/8.5) ‘•’s*'*®] 


b  =  0.38033e°-®®®®2  exp  toossll  +  tanh  [1.49  in  ; 

c  =  in  (9O°/0); 
i  =  cos  [90^  exp(-ac-bc^)]. 


(3.9) 


This  formula  also  has  the  correct  limiting  behavior  when  e->0  or  co,  or  0-^O° 


cr  90°.  An  approximation  formula  for  the  emitted  flux  is  also  available  [Eqs. 
v5j  and  (6)  of  Laframboise  et  al  (1982),  and  Laframboise  and  Kamitsuma 
(1983). 

3.5.  CALCULATION  OF  SECONDARY-ELECTRON  DENSITIES 

Once  the  secondary-electron  escape  fluxes  are  known  (Sec.  3.4),  a  sim.ple, 
inexpensive,  approximate  calculation  of  their  space-charge  density  distribution 
can  be  set  up.  The  proposed  method  is  as  follows;  (1)  ignore  the  gyromotion 
of  the  secondary  electrons  once  they  have  escaped.  Their  motion  then  involves: 
(a)  an  acceleration  along  magnetic  field  lines,  of  amount  -(e/m)E^»B*’/B  (b)  a 

drift  motion  of  velocity  x  across  magnetic  field  lines.  (2)  Integrate 

enough  of  the  trajectories  defined  by  this  motion  (i.e.  their  guiding-center 

trajectories)  to  define  trajectory  tubes  whose  cross-section  at  any  point  can  be 
calculated  with  sufficient  accuracy;  the  method  described  by  Laframboise  et  al 
(1982,  Sec.  7),  can  be  used  to  calculate  the  area  of  a  trajectory  tube  without 
re.f^e.’"ence  to  neighbouring  trajectories.  (3)  Calculate  their  space-charge 
density  n(r^)  at  any  point  by  (a)  ignoring  the  "thermal"  spread  of  their 
velocities  (b)  then  invoking  the  fact  that  their  density  x  their  velocity  [as  given 
by  the  orbit  integration  mentioned  in  (2)],  x  the  cross-sectional  area  A(r  )  of 
the  trajectory  tube  (which  must  be  calculated  in  a  plane  1  the  trajectory)  at 
the  point  r  in  question,  =  a  constant  (whose  value  is  given  by  the  initial 
conditions  at  the  point  on  the  spacecraft  where  the  trajectory  originates)  (c) 
rinding  their  velocity  at  the  point  in  question  by  using  energy  conservation, 
together  with  the  values  of  electric  potential  (p(r)  and  0^  at  that  point  and  the 
emission  point,  and  their  assumed  velocity  v^  at  the  emission  point.  The 


result  is: 


ifr^J  =  n  V  A  /  ■(  Ali^l  V  v  ^  4-  (2e/nn)  [0(r^  1  -  0  ]  I 

O  O  O  I  O  '  /  /  L'T'.  .  Y-qJ  [ 


(3.10) 


where  noVo  is  the  escaping  flux  calculated  in  Sec.  3.4.  At  most  positions,  n(i^ 
will  be  insensitive  to  the  precise  value  assumed  for  Vq^;  assuming  that  Vq  =  the 
cne-sided  therm.ai  speed  (2!<T/7rm)^  will  suffice  for  most  purposes. 
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THETA 

EPS 

15.00 

30 .00 

45. 00 

60 . 00 

75.00 

80 . 00 

85.00 

89.00 

0  .  00 
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.865 

.707 

.500 
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.  174 

.  087 

.017 

.  20 

.991 

.931 

.796 

.585 

.311 

.209 

.105 

.  020 

.50 

.999 

.978 

.892 

.704 

397 

.271 

.  137 

.027 

1  .  00 

.999 

.997 

.971 

.857 

.545 

.384 

.  198 

.  039 

2  00 

.999 

1  .  000 

.999 

.982 

.802 

617 

.342 

.070 

5.00 

.999 

1 .  000 

1  .  000 

1 .000 

.998 

.968 

.723 

.  172 

10  .  00 

.999 

1 .  000 

1  .  000 

.999 

1  .  000 

1 . 000 

.971 

.338 

20 . 00 

.999 

1 .  000 

1 . 000 

.999 

1 .000 

1  .  000 

1.000 

.617 

TABl-E  1 


■^r 


Vol'-^es  cf  the  ratio  i  ~  I/Io  of  ercoping  to  emitted  fl'jx,  fer  varicuo  values 
t?.  ti'o  angle  (in  degrees)  betv/een  the  surface  normal  and  'he  miagnetic-fiold 
direct icr'.,  and  e.  the  nondimensional  repelling  electri  '  field  itrength.  Tr.ese 
‘'VO  '^•jar.tities  apoear  in  the  table  as  TUFITA  and  EPS,  r'^rre';''. r/el v.  These 
results  are  accurate  tc  '.vithin  about  0.5%  or  better:  thus  th.c  differences 
between  .999  and  l.OOG  in  the  Table  are  not  significant.  F  :.r  0  =  0^,  i  =  1 
far  all  values  of  c. 
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'V  ^ 


Total  Ion  Collection  by  a  Unipotentu 
Sphere  in  a  Collisionless  Plasma  Flo> 


Sphere  at  large  negative  potential, 
assuming  that  sheath  thickness 
remains  <  sphere  radius  r.  • 
using  theory  of  Parrot  et  a1  1982/ 
(lower  curve  multiplied  by  / 
1.45)  / 


ion  current 


ion  random 
current 


2.69 


/  \Sphere  at 

/  space  potential, 

/  ”any  Debye  length'^ 

X  . 


2  ' 

/  I 

^  h* —  0.1  eV  protons 


1  1/2^  ,  ,  , 
i  =  — TT  Sj  (ram  current) 


ion  speed  ratio  Sj  =  U/v2kTj/mj 


Figure  2.1  Deoencienco  cf  ion  ‘(j  :o  sp.nere  on  ion  speed  ratio 


-25001 


Ions  O  at  0.1  eV 

Hot  elactrons  at  5000  eV 

Spacacraft  spaad  8  km/aac 


tv: 


Figure  2.2  Deper.der.ee  of  spacecraft  surface  potential  on  hot  electron/total 
ion  a.Tsbient  flux  ratio. 
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Figure  3.1  r.fiect.  cr  surface  orientatic'n  cn  escape  cr  erritted  electrons.  In 
;e',  the  scacecre.^t  surface  is  ce\''per.dicuiar  to  the  magnetic  field  B*",  and 
the  em.itted  electrons,  •Jvnich  exsermnce  an  electric  force  -eE^  directed 
?/vay  from  ‘he  si.'rface,  ail  esceoe.  In  ib),  the  spacecraft  surface  is  nearly 
narallel  to  and  almost  all  •  the  ^':T'.ltted  elec‘rnns  return  to  the  sur 
face,  oven  though  -hey  still  expc-r-tence  can  electric  .^orce  directed  away  frcn 
it.  Note  that  the  component  of  ^  per'pendicular  to  ^  results  only  in  an  E*' 
c  orift  oarailei  to  tne  surtace. 


Figure  3.5.  Ratio  i  =  I/Io  of  escaping  to  emitted  secondary-electron  flux, 
as  a  function  or  the  angle  0  between  the  surface  normal  and  the  magne 
tic  field  direction,  for  varic:un  values  of  the  repelling  electric  field 
strength  param.e^er  r  IF./E'  .t  nv'EkT)^.  The  result  for  e  =  Q  is 
given  bv  i  =  cos  0. 


Nn'e;  [is  preseni 

'9C2.ZB  'X'hcp.  'F.  '  F 


this  program  incorrectly  crecicts  r.cnzerc  electron 
;  h  - 1 ,  i.e.  when  eE*^  has  an  octward  normal 


ment.  out  •.‘■a  oro'ection.  alon.p  nJ  nas  an  inward  normal  ccm.oonent. 


0001  FTN4 

0002  PROGRAM  SCAPE 

0003  C 

0004  C  THIS  PROGRAM  CALCULATES  THE  ESCAPING  FLUX  (NORMALIZED  BY  RANDOM 

0005  C  FLUX)  OF  MAXUELLIAN  ELECTRONS  EMITTED  FROM  A  PLANE  SURFACE  IN 

0006  C  THE  PRESENCE  OF  UNIFORM  ELECTRIC  AND  MAGNETIC  FIELDS.  THE 

0007  C  SURFACE  IS  IN  THE  <X,Y)  PLANE. 

0008  C 

0009  C  DATE  IS  A  iO-CHARACTER  USER-SUPPLIED  DATE  FIELD. 

0010  C  NKSES  IS  THE  NUMBER  OF  CASES  TO  BE  CALCULATED. 

0011  C  NDBUG  IS  A  PARAMETER  GOUERNING  EXTRA  OUTPUT  FOR  DEBUGGING 

0012  C  PURPOSES.  IF  NDBUG  >  0>EXTRA  OUTPUT  OCCURS.  IF  NDBUG  >  1, NDBUG 

0013  C  RATHER  THAN  ONE  ELECTRON  ORBITS  ARE  FOLLOUED . 

0014  C  NALFAjNPSI,  AND  NTH  ARE  THE  NUMBERS  OF  VALUES  OF  ALPHA,  PSI  AND 

0015  C  THETA  (IN  DEGREES)  TO  BE  READ  IN. 

0016  C  NEPS  IS  THE  NUMBER  OF  VALUES  OF  EPS  TO  BE  READ  IN. 

0017  C  NTAU  IS  NUMBER  OF  (DIMENSIONLESS)  TIME  INTERVALS  PER  ELECTRON 

0018  C  GYROPERIOD. 

0019  C  MINT  SELECTS  METHOD  OF  INTEGRATION. 

0020  C  THE  ANGLES  ALPHA  AND  PSI  DEFINE  DIRECTION  OF  THE  ELECTRIC  FIELD 
0021  C  VECTOR  E. 

0022  C  ALPHA  IS  THE  ANGLE  BETWEEN  THE  -E  DIRECTION  AND  THE 

0023  C  SURFACE  NORMAL  (THE  Z  AXIS). 

0024  C  PSI  IS  THE  ANGLE  BETWEEN  THE  PLANE 

0025  C  CONTAINING  THE  2  AXIS  AND  THE  -E  VECTOR,  AND  THE  (X,Z)  PLANE. 

0026  C  THETA  IS  THE  ANCLE  BETWEEN  THE  MAGNETIC  FIELD  VECTOR,  WHICH  IS 

0027  C  ASSUMED  TO  BE  IN  THE  (X,Z>  PLANE,  AND  THE  SURFACE  NORMAL. 

0028  C  EPS  IS  DIMENSIONLESS  ELECTRIC  FIELD  STRENGTH,  DEFINED  AS  THE 

0029  C  POTENTIAL  DIFFERENCE  ACROSS  A  MEAN  EMITTED-ELECTRON  CYRORADIUS, 
0030  C  DIVIDED  BY  THE  EMITTED-ELECTRON  TEMPERATURE  IN  VOLTS. 

0031  C  DVX,OVY,AND  DVZ  ARE  STEPSIZES  FOR  INTEGRATION  OVER  DIMENSIONLESS 

0032  C  VELOCITY.  »«»WARNING»»« : THE  RESULTING  VALUE  OF  NVZ  MUST  BE  NO 

0033  C  LARGER  THAN  THE  DIMENSIONS  OF  VZ , KNECZ , KMIN ,KESC , ZMI N , AND  VZLIM. 

0034  C 

0  035  DIMENSION  VZ ( 1 0 0 ) , KNEGZ ( 1 00 ) ,KMIN( 1 0 0 ) , KESC ( 1 00 ) , ZMIN ( 1 0 0  )  , 

0036  IVZLIM(IOO) 

0037  DIMENSION  X ( 361 ) , Y(361 ) ,SPARL( 361 ) , SPERP ( 361 ) 

0038  DIMENSION  QSI V ( 361 ) ,ETAV( 361 ) , STOR ( 361 ) 

0039  DIMENSION  T AU ( 361 ) ,CSTAU(361 ) ,SNT AU ( 361 > ,Z ( 361 ) , D2 ( 361 ) 

0040  DIMENSION  FLUXdSl) 

0041  DIMENSION  DATE(S) 

0042  DIMENSION  ALPHA ( 91 ) , PSI ( 181 ) , THETA( 181 ) , EPS ( 40 ) 

0043  DIMENSION  ITIME(5) 

0044  C 

0045  PI=3. 14159265 

0046  RTPI=S0RT(PI) 

0047  TOPI=2.0/PI 

0048  TORP=2. 0/RTPI 

0049  READ(7, 10 )DATE , NKSES, NDBUG 

0050  10  FORMAT(5A2,3I5) 

0051  DO  1000  KA3E=1, NKSES 

0052  C 

0053  C  READ  SYSTEM  CLOCK. 

0054  C 

0055  12  ICODE=ll 

0056  CALL  EXEC( ICODE, ITIME) 

0057  FCENS=ITIME( 1 ) 

0058  FSEC=ITIML(2) 

0059  FMIN*ITIME(3) 

0060  FHH-1T1ME(4) 

0061  FDAY=ITIME(5) 

0062  C 
0063 


IF ( NDBUG  GT  0 ) READ ( 7 , 1 20 ) BALPH . BPS  1 , BTHTA , BCP 5 , BVX I N , BVY I N , BVZIN 


0064 
0065 
0066 
0067 
0068 
0069 
0070 
0071 
0072 
0073 
0074 
0075 
0076 
0077 
0078 
0079 
0080 
0081 
0  082 
0083 
0084 
0085 
0086 
0087 
0  088 
0089 
0090 
0091 
0092 
0093 
0094 
0095 
0096 
0097 
0098 
0099 
0100 
0101 
0102 
0103 
0104 
0105 
0106 
0107 
0108 
0109 
0110 
0111 
0112 
0113 
0114 
0115 
0116 
0117 
0118 
0119 
0120 
0121 
0122 
0123 
0124 
0125 
0126 
0127 
0128 
0129 


C 

C 

C 


READ(7,110)NALFA,NPS1,NTH,NEPS,NTAU,MINT 
110  F0RMAT(6I5> 

READ( 7,120 ) (ALPHA(l) ,I=1,NALFA> 
REAO<7,120>(PSI(I) ,I*1,NPSI> 
READ(7,120)(THETA(I),Ii:l,NrH> 
READ(7,120)(EPS(1>,1=1,NEPS> 
READ(7,120>DVX,DVY,DVZ 
120  FORnAT(7E10.3) 

1ESC»12  GIVES  FORK  FEED  ON  H-P  PRINTER. 


IESC»12 

UR ITE ( 6 , 1 25 ) lESC , DATE , K ASE 

125  F0RHAT(A2,1X"FLUX  OF  SECONDARY  ELECTRONS  ESCAPING  FROM  A  SURFACE  I 
IN  UNIFORM  E  AND  B  FIELDS.  ••,5A2,*  ,  CASE^IS) 
URITE(6,130)NALFA,NPSI,NTH,NEPS,NTAU,MINT,DVX,DVY,DVZ 
130  FORMAT</“  NALFA  NPSI  NTH  NEPS  NTAU  MINT"  ,7X“DVX“ ,  7X"DVY"  ,  7X"DVZ  *• 

1  /lX6I5,iP3E10.3) 

URITE<6,i40)(ALPHA<l>,l-l,NALFA> 

URITE(6,i41)(PSI<I),l>l.NPSI> 

URITE<6,142>(THETA<I>,I>1,NTH> 

URITE<6,143) (EPS(I),I«1,NEPS> 

140  FORMAT</-  ALPHA", 1P12E10 . 3) 

141  FORMAT(/"  PSI",1P12E10.3> 

142  FORMAT </-  THETA", lPi2Ei 0 . 3) 

143  FORMAT</"  EP8",iP12E10.3> 

KPBUCbQ 

KUEER«0 

NVX-1NT(9. 0/DVX>>l 

NVY«INT(9.0/DVY)<«-1 

NVZ-INT(4.5/DVZ>4t 

DTAU«2. 0*PI/NTAU 

IF  <  NOBUC . CE . 2 ) NTAU>NDBUG«NT AU 

NTAUP-NTAU^l 

DO  160  ITAU«2,NTAUP 

TAU  < I TAU ) >DT AU»  < I T AU-1 ) 

CSTAU<ITAU)-C08(TAU<ITAU)> 

160  SNTAU<ITAU>-SIN(TAU(ITAU>> 

DO  161  1VZ-1,NVZ 

161  VZ(IVZ)aOVZ«(IVZ-l) 

DO  900  IALFA«1 , NALFA 
ALFA-ALPHA< IALFA>«PI/10B  0 
DO  900  lPSl-l,Nf81 
PS1I«PSI<IPSI)«P 1/180  0 
COPSI>COS<PSII ) 

SIPS1<‘SIN(PSII) 

URITC<6,166>ALPHA(IALFA>  .PSKIPSl) 

166  F0RMAT(//1X"ANCLE8  DEFINING  ELECTRIC  FIELD  DIRECTION  ALPHA  «  " 

A  1PE10.3,"  DECREES,  PSI  »  "IPEIO  3,-  DECREES") 

DO  800  ICPS*1,NEPS 
EPSS>EPS< lEPS) 

EX— EPSS«SIN(ALFA)«COPSl 
EY— EPSS*SIN<ALFA)«SIPSI 
EZ— EPSS«COS(ALFA) 

DO  700  ITH»1,NTH 
FETA-THETA(1TH>«PI/180 . 0 
COSTH>COS(FETA) 

S1NTH*SIN<FETA) 

EQSI>EX*SINTH^CZ«COSTH 
EETA— EXyCOSTH^EZOSINTH 
COA*-EQSl/Pl 
COE=TOPI*EY 
COF-TOPIOEETA 

SUM«U  0 


38 


0130  DO  60U  IOX>l,NUX 

0131  VXIN=-4 . 5+DVX*< IVX-1) 

0132  VXINS=VXIN«eSINTH 

0133  VXINC«VXIN«COSTH 

0134  DO  600  1VY»1,NVY 

0135  VYIN='4.5+DVY»<1VY-1) 

0136  DO  500  IVZ=1,NVZ 

0137  VZIN=VZ<IVZ) 

0138  VQSINaOXINS-»-OZIN«COSTH 

0139  WETIN=-VXINC+yZIN«SINTH 

0140  CODsOQSlN»TORP 

0141  C0C=T0RP*VETIN+T0PI»EY 

0142  COD=TORP*VYlN-TOPI*EETA 

0143  ZMIN<lV2)=>-0.0 

0144  KNEGZ(IOZ)-0 

0145  KHIN(1VZ)=0 

0146  X(1)=0.0 

0147  Y<1)=0.0 

0148  SPARL<l)-0.0 

0149  SPERP<1)=0.0 

0150  Z(l)-0.0 

0151  DQSI=COB 

0152  DETA=COC-COE 

0153  DZ(  1  )»DQSI«CQSTH-«-DETA»SINTH 

0154  D2ZIN=2. 0*COA»COSTH+COD»S1NTH 

0155  IFCVZIN.EQ. 0 . 0 . AND . D2ZIN . LT . 0 . 0 > KNEGZ ( I VZ )»2 

0156  C 

0157  C  SEARCH  ALONG  AN  ORBIT  FOR  A  LOCAL  HINIMUH  IN  Z. 

0158  C 

0159  IF<mNT  GT  .  0)GO  TO  360 

0160  C0G»2. 0«COA«COSTH 

0161  COH=COB*COSTH-COE»SINTH 

0162  COI=COC*SINTH 

0163  COJi^CODKSINTH 

0164  DO  350  lTAUa2,NTAUP 

0165  QSI=>(COA»TAU(ITAU)-^COB)»TAU<ITAU> 

0166  ETA=COC*SNTAU<ITAU)i-COD»(l . 0-CSTAU( ITAU) )-COE«TAU( ITAU> 

0167  X( ITAU)=QSI*SINTH-ETA*COSTH 

0168  Y(ITAU)=COD«SNTAU(ITAU>+COC«(CSTAU(ITAU>-l . 0)+COF*TAU(1TAU> 

0169  SPARL(ITAU)=X(ITAU)»COPSI+Y(ITAU)*SIPSI 

0170  SPERP(ITAU}=-X(ITAU>»SIPSI-fY<ITAU)«COPSI 

0171  Z< ITAU)=aSI»COSTH+ETA»SINTH 

0172  350  DZ< ITAU )=COG*TAU( ITAU) ♦COH+COI»CSTAU( ITAU >+COJ»SNTAU( ITAU) 

0173  CO  TO  380 

0174  C 

0175  C  THE  NEXT  25  STATEMENTS  USE  THE  H-P  "VECTOR  INSTRUCTION  SET" 

0176  C  (ARRAY  PROCESSOR)  TO  REPLACE  LOOP  350,  JUST  ABOVE,  IN  ORDER 

0177  C  TO  SPEED  EXECUTION. 

0178  C 

0179  360  CALL  VSMY ( COA , TAU <2) , 1 ,STOR (2 ) , 1 ,NTAU ) 

0180  CALL  VSAD(C0B,ST0R(2) ,1,ST0R(2) ,1,NTAU) 

0181  CALL  UMPY(ST0R<2),1,TAU<2),1,QSIV(2),1,NTAU) 

0182  C 

0183  CALL  VSSBd  0,CSTAU<2),1,STOR(2),1,NTAU) 

0184  CALL  VSMY<C0D,ST0R(2),1,ST0R(2),1,NTAU) 

0185  CALL  VPIV<C0C,SNTAU<2),1,ST0R(2),1,ST0R<2) ,1,NTAU) 

0186  CALL  VPIV(-C0E,TAU<2),1,ST0R(2) ,1,ETAV(2) ,1,NTAU) 

0187  C 

0188  CALL  VSMY<-C0STH,ETAV(2) , 1 ,ST0R(2) , 1 ,NTAU) 

0189  CALL  VPIV<SINTH,QSIV<2),1,ST0R(2),1,X(2),1,NTAU) 

0190  C 

0191  CALL  VSAD<-1  0,CSTAU<2),1,STOR(2),1,NTAU) 

0192  CALL  VSMY(C0C,ST0R(2),l,ST0R(2),l,NrAU) 

0193  CALL  VPIV<C0D,SNrAU<2) , 1 ,STOR<2) , 1 ,STOR<2) , 1 ,NTAU) 

0194  CALL  UPIV<C0F,TAU<2),1,ST0R<2),1,Y<2),1,NTAU) 

0195  C 


CALL  VSMY <SIPSI,Y<2),1,ST0R<2),1, NTAU ) 

CALL  VPIV<C0PSI,X<2) ,1,ST0R<2) , 1 ,SPARL<2> , 1 ,NTAU) 


0196 
0197 
0198  C 

0199  CALL  VSMY(C0PSI,Y<2),l,ST0R(2),i,NTAU) 

0200  CALL  VPIV(-SIPSI,X<2),1,ST0R(2),1,SPERP(2),1,NTAU> 

0201  C 

0202  CALL  VSMY(SINTH,ETAV<2) ,1,ST0R(2) ,1,NTAU> 

0203  CALL  V'PIV(C0STH,QSlU(2>>l,ST0R(2>,i,Z(2>,l,NTAU) 

0204  C 

0205  COGs2. 0«COA«COSTH 

0206  CALL  USnY(C0C,TAU(2),i»ST0R(2),l,NTAU> 

0207  COH-COB»COSTH-COE«SINTH 

0208  CALL  USAD<C0H,ST0R<2),1,ST0R(2>,1,NTAU) 

0209  CALL  VPlV<C0C#SINTH,CSTAU(2),i,ST0R<2),i,ST0R<2),l,NTAU) 

0210  CALL  VPIV(C0D«SINTH,SNTAU<2>,l,ST0R(2>,i,DZ<2>,i,NTAU) 

0211  C 

0212  380  00  400  ITAU=2,NTAUP 

0213  IF(Z(ITAU)  LT. 0 . 0 )KNEGZ( IOZ>=MAXO (KNECZ< 1UZ> , 1 > 

0214  IF<DZ<ITAU).GE.O.O.AND.DZ(lTAU-i).LT.O.O>GO  TO  443 

0215  400  CONTINUE 

0216  1F(KNECZ(IUZ)  GT.OIGO  TO  402 

0217  C 

0218  C  NO  NEGATIVE  VALUES  OF  Z  OR  HINIHA  IN  Z  HAVE  BEEN  FOUND 
0219  C 

0220  401  KESCdVZ)-! 

0221  GO  TO  470 

0222  C 

0223  C  NEGATIVE  VALUES  OF  Z  HAVE  BEEN  FOUND,  BUT  NO  HINIHUM  IN  Z  PRINT 

0224  C  ORBIT  FOR  EXAMINATION. 

0225  C 

0226  402  KESC<IVZ>*=0 

0227  1F<KUEER .CE.3)GO  TO  470 

0228  URITE(6,409) 

0229  409  FORMAT</"  ORBIT  DETAILS") 

0230  URXTE(6,410) 

0231  410  FORMAT <6XSHALPHA,7X3HPSI,SXSHTHETA,6X4HEPSS,6X4HVXIN,6X4HVYIN, 

0232  A  6X4HVZIN) 

0233  UR1TE(6,411)ALPHA(IALFA>,PS1(IPSI> ,THETA(1TH> ,EPSS,VXIN, VYIN, VZIN 

0234  411  FORMAT(1X1P7E10.3) 

0235  URITE(6,412)(I,X(I>,Y<1) ,Z(I>,DZ(I>,I-1,NTAUP> 

0236  412  F0RMAT</1X3<5X"I",5X"X<I)“,5X"Y(I)",5X"Z<I)",4X"DZ(I> ")/ 

0237  A  <1X3(I6,4F9  3)) ) 

0238  KUEER=KUECR>1 

0239  GO  TO  470 

0240  C 

0241  C  A  MINIMUM  IN  Z  HAS  BEEN  FOUND.  DETERMINE  ITS  Z  VALUE  ZMlNdVZ) 

0242  C 

0243  443  KMINdVZ)  =  l 

0244  FRACT  =  DZdTAU-l  )/<DZdTAU-l)-DZdTAU)  ) 

0245  TAMIN=TAUdTAU-l  )+FRACT*DTAU 

0246  aSI=(COA»TAMIN>COB)«TAMIN 

0247  ETA=COC*SIN<TAMIN)+COD«<l ,0-COS{TAM1N) )-COE«TAMlN 

0248  ZMIN< IVZ)=QSI»COSTH+ETA*SINTH 

0249  IF<ZMINdVZ)  LT  .  0 . 0  )KMINd  VZ  )  »2 

0250  KESCdVZ)=0 

0251  IF<ZMINdVZ)  GE  .  0  0>KESCdVZ>»l 

0252  C 

0253  470  IF(NDBUG  EQ.OGO  TO  500 

0254  IF(BALPH  NE  ALPHA < I ALFA >) GO  TO  500 

0255  IF<BPSI  NE  PSIdPSX))GO  TO  500 

0256  IF<BTHTA.NE.THETAdTH))GO  TO  500 

0257  IF(BEPS.NE  EPSS)GO  TO  500 

0258  IF<ABS<VXIN-BVXIN)  GT  O.OOODGO  TO  500 

0259  IF(ABS<VYIN-BVYIN)  GT  O.OOODGO  TO  500 

0260  IF<ABS<VZIN-BVZIN)  GT.O  OOODGO  TO  500 

0261  WRITE(6,409> 


0263  URITe(6,410) 

0263  WRITE(6,411 )ALPHA< lALFA) ,PSI ( IPSI ), THETA < ITH >, EPSS , VXIN , VYIN , V2IN 

0264  URITE(6,412)<I,X(1>,Y(1>,Z(1),DZ(1),I=1,NTAUP) 

026S  500  CONTINUE 

0266  C 

0267  C  KESC(IUZ)=i  OR  0  DEPENDING  ON  WHETHER  THE  lUZ'TH  ORBIT  DID  OR 

0268  C  DID  NOT  ESCAPE.  KHIN(IUZ)  >  0  OR  >  0  DEPENDING  ON  WHETHER  A  LOCAL 

0269  C  HINIMUH  IN  Z  WAS  OR  WAS  NOT  FOUND.  KLIM  IS  THE  FOUND  NUHBER  OF 

0270  C  STARTING  OR  END  POINT  VALUES  VZL1M( i ,2, . . . ,KLIN)  FOR  INTEGRATION 

0271  C  IN  UZIN  (THE  INITIAL  Z  VELOCITY),  FOR  ESCAPING  ORBITS. 

0272  C 

0273  KLIM=0 

0274  IF<KESC(1)  EQ. 0)CO  TO  501 

0275  KHM=1 

0276  VZLIM(1)=0.0 

0277  501  DO  510  IVZ>2,NVZ 

0278  IF<KESC(IVZ)  .E(J.KESC<IVZ-l>)CO  TO  510 

0279  KLIM=XLIM+1 

0280  IF(KMIN(IVZ-1) .EQ.O.OR.KHIN<IVZ>.EQ.O>CO  TO  505 

0281  VZLIH(KLIM)aVZ<IVZ-l)+ZHIN<IVZ-l>/(ZMIN(IVZ-l)-ZMIN<IVZ) >* 

0282  A  <VZ<IVZ)-VZ<IVZ-1)  ) 

0283  CO  TO  510 

0284  C 

0285  C  AN  INTERVAL  IN  VZIN  HAS  BEEN  FOUND  CONTAINING  A  CHANGE  BETWEEN 

0286  C  ESCAPE  AND  NO  ESCAPE,  BUT  A  LOCAL  HINIMUH  IN  Z  IS  NOT  FOUND 

0287  C  FOR  ONE  OR  BOTH  ENDS  OF  THIS  INTERVAL.  PRINT  ORBIT  PARAMETERS  FOR 

0288  C  EXAMINATION. 

0289  C 

0290  505  VZLIM(KLIM)=0.5»<VZ(IVZ-i)+VZ<IVZ>) 

0291  KPBUG=KPBUG+1 

0292  IF<KPBUG.GE.6.0R.NDBUG.EQ.O>CO  TO  510 

0293  WRITE(6,409) 

0294  WRITE<6,410) 

0295  WR ITE (6,411)  ALPHAdALFA)  ,PSI(  IPSI ), THETA ( ITH)  ,EPSS,VXIN, 

0296  1  VYIN,VZ(IVZ) 

0297  C 

0298  510  CONTINUE 

0299  C 

0300  SIGN=i.O 

0301  VZSUM=0.0 

0302  IF(KLIM.Eq. 0)GO  TO  530 

0303  DO  520  ILIM=1,KLIM 

0304  VZL=VZLIM< ILIM) 

0305  VZSUM»VZSUM  +  SICN»EXP ( -V2L*VZL  > 

0306  520  SIGN=-SICN 

0307  C 

0308  530  IF(NDBUG.EQ  0>GO  TO  600 

0309  IF(BALPH  NE. ALPHA( lALFA)  OR .BPSl ,NE  PSK IPSI) )GO  TO  600 

0310  IF(BTHTA.NE.THETA(ITH) .OR.BEPS.NE.EPSS)CO  TO  600 

0311  IF(ABS<VXIN-BVXIN)  GT. 0001, OR. ABS(VYIN-BVYIN> .GT. . 0001)00  TO  600 

0312  WRITE(6,540 ) ( I , KESC ( I ) , KMIN( I ) , ZMIN ( I > , 1=1 ,NVZ) 

0313  540  F0RMAT(/1X4( "  I  KESC(I)  KMIN(I)  ZMIN<I)")/ 

0314  A  ( 1X4( 15,218, 1P310 . 3) ) ) 

0315  IF(KLIM.GT . 0 ) WR I TE ( 6 , 550 ) ( VZLIM( I ) , 1 =l ,KL IM ) 

0316  550  FORMAT(/("  VZLIM  "iP12E10.3)) 

0317  URITE(6,S60)VZSUM 

0318  560  FORMAK/"  VZSUM  =  "IPEIO  3) 


SUM=SUM^VZSUM*EXP(-VX1N*VXIN-VYIN*VYIN) 


700  FLUX(  ITH)=SUM*0VX4DVY/PI 
IFdEPS  GT  DUO  TO  771 
WRITL(6,768) ( THETA (ITH) ,ITH=1 ,NTH) 

768  FORMAr</9X"THETA",14Fa4/(20X14F84)) 
URnE(6,770) 

770  FORMAT  (■•  EPS  "/) 


EPS 


0328  771  WRITE(6,772)EPSS,<FLUX{ITH>,ITH*1,NTH) 

0329  772  FORrtAT<lXlPE10.3,3X0P14F8.5/(20X0P14F8.5)) 

0330  C 

0331  800  CONTINUE 

0332  900  CONTINUE 

0333  C 

0334  C  READ  SYSTEM  tLOCK  AGAIN  AND  PRINT  ELAPSED  TIME. 

033S  C 

0336  CALL  EXEC( ICODE, ITIME) 

0337  FC2=ITIME<1) 

0338  FS2«=ITIME<2) 

0339  FM2=IT1ME<3) 

0340  FH2=ITIME<4) 

0341  FD2=ITIME<S) 

0342  TIMIN»<24 . 0« < FD2-FDAY >-i-FH2-FHR  >»60 . 0+FM2-FmN-KFS2-FSEC^0 . 01*(FC 

0343  A  FCENS) )/60 . 0 

0344  URITE<6,910)KPBUC,TIMIN 

034S  910  F0RMAT</63XIS,-  INTEGRATION  LIMIT  FIXUPS.  EXECUTION  TIME  ", F8.2, 

0346  A  “  MINUTES") 

0347  1000  CONTINUE 

0348  END 

0349  END« 


